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Abstract. We report a measurement of the CP violation in S" — > J/if/K^ decays. We perform a time-dependent 
analysis of the decays reconstructed in 35pi? ' of LHCb data that was taken in 2010. We measure the CP asym- 
metry parameter Sj/^go = 0.53*^-^ (stat) ±0.05 (syst), which is coimected to the CKM angle jS. 



1 Introduction 

The decay B" J/tfrKg is well known as the gold-plated 
mode for the study of CP violation in the Z?" meson decays 
to a CP eigenstate common to both Z?" and allowing for 
interference through oscillation. Therefore, measurements 
of the decay and have a good sensitivity to the CKM 
angle sin{2p), which is connected to the parameters S j/^jq 
andCy/^xo: 

= S ji^Kl siniAnidt) - Cj^^jq co%{Amdt) , (1) 

where zfm^ = m^o -m^o is the mass difference between 
the mass eigenstates. We have neglected the decay width 
difference between the mass eigenstates of the system. 
The coimection of the S ji^^o and Cjj^j^ parameters to the 
CKM angle is 

S ji^Kl - >/l - Cji^go^ sin . (2) 

In the Standard Model, both CP violation in mixing and 
direct CP violation are negligible in — > ccs decays. As a 
consequence the cosine term vanishes, implying 

Sji^Kl^^^W- (3) 

During the last decade the B-factories BABAR and Belle 
reached outstanding precision in the measurement of S jj^jm . 
The most recent BABAR measurement reports S j/^go = 
0.663 + 0.039 (stat) ± 0.012 (syst) fl]. The most recent 
Belle result is S j/^jq = 0.642 + 0.031 (stat) ±0.017 (syst) 
[2]. Currently, the world average [3] is sin2p = 0.673 ± 
0.023. 

The measurement presented in this note [4] with the 
first LHCb data is an important step to demonstrate the po- 
tential of the LHCb experiment in this topic and it demon- 
strates that the flavour tagging algorithms are under con- 
trol. We measure S j/^j^a under the assumption that C j/^fq = 
0, and we quote the resulting value of C j/^j^o if the Stan- 
dard Model constraint is relaxed. 
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2 Analysis Strategy 

The results presented are based on data collected with the 
LHCb detector at the LHC collider at a center-of-mass en- 
ergy of V* = 7 TeV. Details on the LHCb experiment can 
be found in elsewhere [5]. The data set analyzed has an 
integrated luminosity of « 35pb"^ We use Monte Carlo 
(MC) simulated samples that are based on the Pythia 6.4 
generator [6]. The EvtGen package [7] was used to gener- 
ate hadron decays and the Geant4 package [8] for detector 
simulation. 

We reconstruct J/if/ candidates in the decay mode J/iJ/ 
IJ-^jj.', from pairs of opposite sign tracks that have a trans- 
verse momentum of pj > 500MeV/c each and particle 
identification signatures consistent with those of muons. 
The invariant mass of the pair must be compatible with the 
known Jlij/ mass. K^^ candidates are reconstructed through 
their decay into n^n' with take pairs of oppositely charged 
tracks. An additional requirement of LIctl > 5 is made on 
the candidate, where L is the distance between the K'^ 
decay vertex and the B decay vertex and ctl is the uncer- 
tainty of L. We constrain the invariant masses of the recon- 
structed y/i/' and candidates to their known masses. 

The initial B flavour is determined by the combination 
of various tagging algorithms. These either determine the 
flavour of the non-signal b hadron produced in the event 
{opposite side, OS), or they search for an additional pion 
accompanying the signal Z?" or Ir' (same side, ^Sn). There 
are four tagging algorithms that use the charge of the lep- 
ton (ju, e) from semileptonic B decays, or that of the kaon 
from the b c ^ s decay chain, or the charge of the in- 
clusive secondary vertex reconstructed from b decay prod- 
ucts. All of these algorithms have an intrinsic mistag rate, 
due to picking up tracks from the underlying event, or due 
to flavour oscillations of neutral tag B mesons. For each 
signal Z?" candidate the tagging algorithms also predict the 
mistag probability o). For this various kinematic variables 
such as momenta and angles of the tagging particles are 
combined into neural networks. The neural networks are 
trained on MC simulated events. 

The flavour asymmetry that is accessible in Z?" — > Jjij/K^ 
decays directly depends on the dilution D due to the mistag 
probability for each signal candidate, {d^"^ - ^ 2,(1 - 
2a),)^. Its statistical precision is proportional to the inverse 
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square root of the effective tagging efficiency eeff> 

eeif = etag , (4) 

where Stag is the probabiHty that a tagging decision is found. 
The tagging algorithms are optimized for highest Seff on 
data, using the self-tagging decays — > Jjij/K^ and 
D*'ji^v. The estimated mistag probability is calibrated on 
these same channels. The effective tagging efficiency is 
measured to be £eff = (2.82 ± 0.87)%. 

To extract the S j/^go parameter, we perform a simul- 
taneous, extended unbirmed maximum likelihood fit to the 
proper time and the invariant mass distributions: 

S 1=1 

We minimize - In X to find optimal values for the fit 
parameters A. 

The fit is simultaneous in four subsamples s. These 
subsamples are defined by whether the candidates were 
triggered by a lifetime unbiased ("U") or a lifetime biased 
("B") decision; and by whether or not a tagging decision is 
available ("u" for untagged, "t" for tagged). Each subsam- 
ple contains events, and « = 2* 

We consider four observables: the reconstructed mass 
m of the candidate (5.15 GeV/c^ <m< SAGeV/c^), its 
proper time ? (- 1 ps < / < 4 ps), the flavour tag decision 
d, and the combined per-event mistag prediction o). The 
flavour tag d can take the discrete values J = 1 if tagged as 
initial and J = -1 if tagged as initial B". 

The probabihty density functions (p.d.f.s) consist 
of three components, signal (S), prompt background (P), 
and long lived background (L). The mass p.d.f. of the sig- 
nal component consists of a single Gaussian. We assume 
both background components have similar mass distribu- 
tions, and we use the same parameterization for both. It 
is modeled as an exponential with a single shape param- 
eter. The proper time p.d.f. of the signal component can 
be written as Psit,d,a)) = Psit,d\a)) ■ fsioh- The first 
term is a conditional p.d.f. as it depends on the value of 
o), the second term describes the distribution of w. The 
background parameterization of the proper time factorizes, 
feit, d, CO) = feit, d) ■ faioS), where B = P,L. 

In the fit to the B" J/if/K^ channel we fix the mix- 
ing frequency Anid to its nominal value of Anid - (0.507 ± 
0.005) • 1012 hs-^ [9], and Cj/^,q = 0. In total, there are 

27 floating parameters: the CP parameter S j/^jq, the B" 

hfetime t, the mass Wq , twelve event yields, four pa- 
rameters of the long-lived proper time background, five pa- 
rameters of the time resolution, the mass signal resolution 
cr^, and two parameters of the mass background shape. We 
have checked the fit implementation on a large sample of 
MC generated signal events, where we find good agree- 
ment with the generated values. 

The measurement was performed using a "blind" anal- 
ysis technique to minimize unconscious experimenter bias, 
the parameter of interest was encrypted in the likelihood 
fit. Only after the full analysis strategy was developed and 
proved to be stable, the encryption was removed, unbhnd- 
ing the true result. 



Table 1. Fit result of the nominal fit to the full J/i//K^ data 
sample. 



Parameter 


Unit 


Fitted Value 


nts 

0-S,m 
T 


MeV/c^ 
MeV/c^ 
ps 


5278.13 ±0.29 
8.82 ± 0.24 
1.517 ±0.046 


al 
Jlj 

L,l 

/r,2 
CK3 


(MeV/c^r' 
(MeV/c^)-' 

ps 
ps 
ps 

ps 
ps 
ps 


-8.71 ±3.8- lO^"* 
-5.86 ± 0.87 • 10"* 
0.836 ± 0.054 
0.221 ± 0.036 
1.04 ±0.24 
0.482 ± 0.029 
0.500 ± 0.019 
0.477 ± 0.017 
0.02522 ± 0.00066 
0.0685 ± 0.0016 
0.293 ± 0.019 



2.1 Results 

The result of the maximum likelihood fit to the full data 
sample is summarized in Tables 1 and 2. The mass and 
proper time distributions and the fit projections are shown 
in Figure 1 . Figure 2 shows the resulting time dependent 
raw asymmetry, which contains all fit components. The 
asymmetry in the lowest proper time bins is therefore dom- 
inated by the backgrounds, whereas the measured asym- 
metry in the high proper time bins is dominated by signal 
events. The measured value of 5 j^^j^ is 

^ j/^K° = 0.53;o;29 , (6) 

where the error is statistical only. We find the global cor- 
relation coefficient of S j/^/q to be piS j/^iq) - 0.016. We 
also perform the nominal fit without the Standard Model 
constraint Cjj^^ = 0. In this case, we find 

Cj,^^^ = 0.28 ± 0.32 , S = 0.38 ± 0.35 , (7) 

again quoting statistical errors only. The correlation coef- 
ficient between the parameters is p(S j/^[^o, Cj/^jq) = 0.53. 
Their correlations to other parameters are neghgible. 

Table 2. Fitted event yields in the full ^ J/iJ/K^ data sample. 



Sample 


Parameter 


Fitted Value 






221 ± 17 


U,t 




3218 ±62 






309 + 33 






59.8 + 8.7 


B,t 




164 ± 14 






102 ± 12 






767 ± 32 


U,u 


21134 ± 161 






807 ± 79 






279 ± 18 


B,u 


747 ± 30 






339 ± 23 
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\f^ = 7TeV L = 35pB' 
3/\\l K 




Statistical significance of a non-zero CP violation from the 
likelihood ratio of a test fit, in which we fix S j/^^ - 0, to 
be 



S = ^/2]n(X^^^^^JX~n) = 1.8 



(9) 



This is the first CP violation result in the golden chan- 
nel 5°^ in LHCb. 
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Fig. 1. Reconstructed mass (left) and proper time (right) distribu- 
tions of B" — > J/tf/K'^ candidates. Overlayed are projections of the 
component p.d.f.s used in the fit: full p.d.f. (solid black), signal 
(dashed blue), prompt background (dash-dotted red), long lived 
background (dotted orange). 
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Fig. 2. Time dependent raw CP asymmetry in B" — » J/if/K^. The 
solid curve is the full p.d.f. (signal and background) overlayed 
onto the data points. The green band corresponds to the one stan- 
dard deviation statistical error. 



3 Conclusion 

The final result on the CP violation parameter 5 y^^^^ is 

Sj/^j,n = 0.53!°:g(stat) + 0.05(syst) . (8) 

This result is compatible with the World Average, and 
dominated by the statistical uncertainty. We calculate the 



